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The chromoelectric and chromomagnetic fields, created by a static gluon-quark-antiquark sys- 
tem, are computed in quenched SU(3) lattice QCD, in a 24 3 x 48 lattice at (3 = 6.2 and 
a = 0.07261(85) fm. We compute the hybrid Wilson Loop with two spatial geometries, one with a 
U shape and another with an L shape. The particular cases of the two gluon glueball and quark- 
antiquark are also studied, and the Casimir scaling is investigated in a microscopic perspective. 
This microscopic study of the colour fields is relevant to understand the structure of hadrons, in 
particular of the hybrid excitation of mesons. This also contributes to understand confinement with 
flux tubes and to discriminate between the models of fundamental versus adjoint confining strings, 
analogous to type-II and type-I superconductivity. 
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Here we present the first Lattice QCD study of the 
chromoelectric and chromomagnetic fields, created by 
a static gluon-quark-antiquark system. Although the 
colour fields have been extensively studied for the quark- 
antiquark, [1-4], for three quarks [4-8], for the hybrid 
only the static potential has been studied so far [9, 10]. 

The hybrid static potential is also relevant to under- 
stand the nature of confinement and of Casimir Scal- 
ing, since with the hybrid potential we can interpo- 
late between the gluon-gluon interaction and the quark- 
antiquark interactions which are particular cases of the 
hybrid static potential. The first study of the static 
gluon-gluon interaction was performed by Michael [11, 
12], and Bali [13] extended this study to other SU(3) 
representations, leading to the Casimir Scaling picture. 
Bicudo et al. [10] and Cardoso et al. [9] studied the 
static gluon-quark-antiquark potential and showed that 
when the segments gluon-quark and gluon-antiquark are 
perpendicular, the potential V is compatible with the 
confinement realized with a pair of fundamental strings, 
one linking the gluon to the quark and the other link- 
ing the same gluon to the antiquark. For parallel and 
superposed segments, however, the total string tension 
becomes larger and is in agreement with the Casimir Scal- 
ing measured by Bali [13]. Bicudo, Cardoso and Oliveira 
established an analogy between the static potential and 
a type-II superconductor for the confinement in QCD, 
illustrated in Fig. 1, with repulsion of the fundamental 
strings and with the string tension of the first topolog- 
ical excitation of the string (the adjoint string) larger 
than the double of the fundamental string tension. In 
type-I superconductor the fundamental strings would be 
attracted and would fuse into an adjoint string. With 
the computation of the flux tubes we can further under- 
stand, microscopically the Casimir Scaling. For instance 
Semay [14] presented a model for Casimir Scaling, based 
on a shape of flux tubes independent of the colour SU(3) 
representation, and we can test it. 

In this paper, we investigate the chromoelectric and 
chromomagnetic fields, and the resulting lagrangian and 
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Figure 1: String attraction and fusion, and string 
repulsion, respectively in type I and II superconductors. 



energy density distributions around a static gluon-quark- 
antiquark system in quenched SU(3) lattice QCD. In sec- 
tion II, we introduce the lattice QCD formulation. We 
briefly review the Wilson loop for this system, which was 
used in Bicudo et al. [10] and Cardoso et al. [9], and show 
how we compute the colour fields and the lagrangian and 
energy density distribution. In section III, the numerical 
results are shown, including several density plots of the 
chromo flux tubes, and longitudinal plots of the chromo 
field profiles. Finally, we present the conclusion in section 
IV. 



II. THE WILSON LOOPS AND COLOUR 
FIELDS 

In principle, any Wilson loop with a geometry similar 
to that represented in Fig. 2a, describing correctly the 
quantum numbers of the hybrid, is appropriate, although 
the signal to noise ratio may depend on the choice of 
the Wilson loop. A correct Wilson loop must include an 
SU(3) octet (the gluon), an SU(3) triplet (the quark) and 
an SU(3) antitriplet (the antiquark), as well as the con- 
nection between the three links of the gluon, the quark 
and the antiquark. 

We construct the gluon-quark-antiquark Wilson loop 
from the two-color-octet meson operator, 

0{x) = \ [q(x)X a r iq (x)j [q(x)X a r 2 q(x)} (1) 

where A a are the Gellmann SU(3) colour matrices, and 
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Figure 2: (a) Wilson loop for the gqq and equivalent 
position of the static antiquark, gluon, and quark, (b) 
Simple Wilson loops that make the gqq Wilson loop. 



where I\ are spinor matrices. Using the lattice links to 
comply with gauge invariance, the second operator in Eq. 
(1) can be made nonlocal to separate the quark and the 
antiquark from the gluon, 



O(x) 



q(x)X a T iq (x) 
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q(x + x 2 jl 2 ) 
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The contraction of the quark field operators, assuming 
that all quarks are of different nature, gives rise to the 
gluon operator, 

1 
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Using the Fiertz relation 
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we can prove that 
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(a) U shape geometry. 



(b) L shape geometry. 



Figure 3: gluon-quark- antiquark geometries, U and L 
shapes. 



where W\ , W 2 and W3 are the simple Wilson loops shown 
in Fig. 2b. Importantly for the study of the Casimir Scal- 
ing, when n = 0, W\ = 3 and W 2 = W3, the operator 
reduces to the mesonic Wilson loop and when /1 = v and 
ri = r 2 = r, W 2 = W\ and W3 = 3, W gqq - reduces to 
W gq q(r,r,t) = \W(r, t)\ 2 — 1, that is the Wilson loop in 
the adjoint representation used to compute the potential 
between two static gluons. 

In order to improve the signal to noise ratio of the 
Wilson loop, the links are replaced by "fat links", 

M*)-^(3)^(M*) 

+wY,U v (s)U yL (s + v)Ut(s + ii)Y (6) 

We use w = 0.2 and iterate this procedure 25 times in 
the spatial direction. 

We obtain the chromoelectric and chromomagnetic 
fields on the lattice, by using, 



(Ef) = (P 0i ) 



(WP 0t ) 
(W) 



and, 



/Q2V _ (WP jk ) 



(7) 



(8) 



where the jk indices of the plaquette complement the 
index i of the magnetic field, and where the plaquette is 
given by 

iV ( s ) = 1 " ^ ReTr P^s)U v (s + fi)Ul(s + v)Ut(s)] . 

(9) 

The energy (H) and lagrangian (£) densities are given 
by 



H 



((E 2 ) + (B 2 )) 



(10) 



(11) 



Notice that we only apply the smearing technique to the 
Wilson loop. 
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Figure 4: Results for the static two gluon glueball. The energy density plot, (c), have a flux tube, responsible for the 
string tension, but we choose to put all the colour scales on the same value in order to be able to make a comparison 
with the values of the different fields, thus the flux tube energy is less visible. The top value of the colour scale is the 
maximum value of the field. The results are in lattice spacing units. 
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(d) Lagrangian Density 



Figure 5: Results for the static quark-antiquark system. The energy density plot, (c), have a flux tube, responsible 
for the string tension, but we choose to put all the colour scales on the same value in order to be able to make a 
comparison with the values of the different fields, thus the flux tube energy is less visible. The top value of the 
colour scale is the maximum value of the field. The results are in lattice spacing units. 



III. RESULTS 

Here we present the results of our simulations with 
266 SU(3) configurations in a 24 3 x 48, (3 = 6.2 lattice, 
generated with the version 6 of the MILC code [15], via 
a combination of Cabbibo-Mariani and overrelaxed up- 
dates. The results are presented in lattice spacing units, 
defined in Eq. (3). 

In this work two geometries for the hybrid system, 
gluon-quark-antiquark, are investigated: a U shape and 
a L shape geometry, both defined in Fig. 3. 

In the U shape geometry, we only change the distance 
between quark and antiquark, d = 0, 2, 4, 6, and we fix 
the distance between gluon and quark-antiquark at I = 8. 
When the quark and the antiquark are superposed, d = 0, 
the system corresponds to a two gluon glueball, Fig. 4. 

In the L shape geometry, we fix the distance between 
the gluon and quark at T2 = 8 and the distance between 
the gluon and antiquark is changed, T\ = 0,2,4,6,8. 
When the gluon and the antiquark are superposed, T\ — 



0, the system is equivalent to a meson. The results for 
the meson system are presented in Fig. 5. 

A. Flux Tube and Casimir Scaling 

First we discuss the results for the two degenerate 
cases, in which the system colapses into a two body sys- 
tem - the meson (L geometry with T\ = 0) and the two 
gluon glueball (U geometry with d = 0). In the meson 
case we confirm the results obtained in previous works 
(for example [2]). Not only in the meson case, but also 
in general, we have (Ed) > (E\) > \ (B]_) \ > \(Bp\. 

We also observe that (E 2 ) > and (B 2 ) < 0, for all 
the studied geometries. Since the absolute value of the 
chromoelectric field dominates over the absolute value 
of the chromomagnetic field, there is a cancelation in en- 
ergy density, Eq. (10), and an enhancement in lagrangian 
density, Eq. (11). 

We measure the quotient between the energy densities 
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Figure 6: Results for the glueball (d = and Z = 8, U 
geometry) energy density over the meson (n = and 
T2 = 8, L geometry) energy density for x=0. Casimir 
scaling, were the Casimir invariant A$ • Xj produces a 
factor of 9/4 (broken line). 



of the meson system and of the glueball system, in the 
mediatrix plane between the two particles (x = 0). The 
results are shown in Fig. 6. As can be seen, these re- 
sults are consistent with Casimir scaling, with a factor 
of 9/4 between the energy density in the glueball and 
in the meson. This corresponds to the formation of an 
adjoint string between the two gluons. The results are 
compatible with an identical shape of the two flux tubes, 
but with a different density, and in this sense this agrees 
with the simple picture for the Casimir Scaling of Semay, 
[14]. 

The results for the fields in the case of the two gluon- 
glueball are given in Fig. 4 and the meson in Fig. 5. 

B. L Geometry 

The squared field components in the L geometry with 
7*1 = T2 = 8 are shown in Fig. 7. In this figure, we can 
see that (E 2 ) is greater is the x axis and (Ey) is greater 
in the y-axis, on the other hand the chromomagnetic field 
components exhibit the reciprocal behaviour - | (B%) | is 
greater in the y axis and | (By) | is greater on the x axis. 
This result is consistent with having two essentially in- 
dependent fundamental strings, since this was the result 
obtained for one fundamental string - the longitudinal 
component is the dominant one in the chromo-electric 
field and the transversal component is dominant in the 
chromo-magnetic field. 

In Fig. 9e-h , we show the distribution of the la- 
grangian density, in the L geometry, with 7*2 = 8, 
fixed, and for different ri, where r\ is the distance 
between gluon-antiquark and r<i the distance between 
gluon-quark. The variation of the lagrangian density 
with ri can also be seen in Fig. 8c and in Fig. 8d, in the 
x and y axis (Fig. 3), where is the anti-quark and the 
quark. Notice that the result in the y axis is essentially 
the same, when we move the antiquark in the x axis, ex- 
cept for the case of r\ =0, where the system collapses in 



a meson. But, even in this case, the flux tube near the 
quark is almost the same. 

In the y axis, we can see the presence of a flux tube 
between the gluon and the quark. As can be seen for T2 — 
8, the lagrangian density tends to a constant in the center 
of the tube and remains practically unchanged when the 
antiquark and the gluon are far apart. This last result 
is consistent with the existence of a confining potential 
V gq — > err between the gluon and the (anti) quark. 

Our results indicate that in this geometry the system 
is well described by two independent fundamental strings 
as was stated in [10] and [9]. 

C. U Geometry 

We show the results for the chromoelectric and chro- 
momagnetic fields in the U geometry at distances I = 8 
and d = 6 in Fig. 10. The results are consistent with 
the ones for the L shape geometry. The longitudinal 
component of the chromoelectric field is the dominant 
component. This is the y component of the chromoelec- 
tric, and this is expected since the flux tube is essentially 
aligned in this direction. In the same way (B%) and (B 2 Z ) 
are seen to be dominant with relation to (By), which is 
consistent with the fact that the transversal component 
of the magnetic is the larger one. 

In Fig. 9a-d we can see the evolution of the lagrangian 
density, as a function of the the quark- antiquark distance, 
d, for fixed 1 = 8. For d = 0, we are in the glueball 
case and we thus have an adjoint string linking the two 
gluons. For d = 2 we can see the stretching of the tube 
in the x direction. For d = 4, corresponding to y ~ 5, 
we already see the string splitting in two fundamental 
strings. In d=6, the separation of the two fundamental 
strings is clear, they only join at the gluon position. The 
transition point between the two regimes - one adjoint 
string and two clearly splitted fundamental ones - occurs 
between d = 2 and d = 4, for I = 8. This transition 
point occurs for an angle between the two fundamental 
strings of 0.37 ± 0.12 rad, and this is relevant for the 
quark and gluon constituent models. In Fig. 8a, we can 
see the stretching and partial splitting of the flux tube 
in the equatorial plane (y = 4) between the quark and 
the antiquark, and in 8b we see the results for the y axis, 
where the gluon is located (at y = 0), as well and the 
centroid of the qq subsystem (at y = 8). 

IV. CONCLUSIONS 

We present the first study the chromoelectric and chro- 
momagnetic fields produced by a static quark-gluon- 
antiquark system in a pure gauge SU(3) QCD lattice. 

We report the cases of a simple meson and of a two 
gluon glueball, which correspond to two different degen- 
erate cases of a hybrid meson system. We verify the qual- 
itative results for the squared components of the colour 
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Figure 7: Chromoelectric and chromomagnetic components and energy and lagrangian densities in the L shape 
geometry for ri = 8 and r2 = 8. We use different colour scales to have a better view of the flux tube and the top 
value of the scale is the maximum value of the field. The results are in lattice spacing units. 



fields that were obtained by other authors for the quark- 
ant iquark system. Namely, we find that the chromoelec- 
tric field is dominant over the chromomagnetic and that 
the longitudinal components of the chromoelectric field, 
as well as the transversal component of the chromomag- 
netic, are dominant over the other components of the 
respective fields. We find a similar behaviour in the two 
gluon system. We also verify that the results for the 
two degenerate systems are related, with the energy den- 



sity in the glueball flux tube being compatible with 9/4 
times the energy density in the meson flux tube. This 
is in agreement with the Casimir scaling factor between 
the glueball and the meson, obtained by Bali [13]. 

We also study two geometries for the hybrid meson 
system. We study a L shaped geometry, with the gluon 
on the origin, the quark on the y axis and the antiquark 
on the x axis. In this case we verify the dominance of 
the longitudinal component in the chromoelectric field 
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(a) U geometry at y = 4. 



6 8 



r1 =0 and r2=8 
r1 =2 and r2=8 
r1 =4 and r2=8 
r1 =6 and r2=8 
r1 =8 and r2=8 




0.0e+00 



(c) L geometry along the segment gluon-antiquark. 
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(b) U geometry at x = 0. 
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(d) L geometry along the segment gluon-quark. 



Figure 8: Results for the lagrangian density. The lines were drawn for convenience and therefore do not represent 
results from any kind of interpolation. The results are in lattice spacing units. 



and of the transversal component in the chromomagnetic 
field in the two flux tubes coming from the gluon. We 
also concluded that this two flux tubes are, mainly, two 
independent fundamental strings, which agrees with the 
results for the potential obtained by [10] and [9]. We also 
study a U shaped geometry, which allow us to see the 
transition between the two regimes of confinement, with 
one adjoint and with two splitted fundamental strings. 

Whether the Casimir scaling is due to a repulsive su- 
perposition the two fundamental strings or to the actual 
existence of an adjoint string, we cannot yet distinguish 
in the present study. We conjecture that both these two 
pictures are essentially equivalent in the gluon-gluon sys- 



tem. But it appears that for angles between the gluon- 
quark and gluon-antiquark segments larger than 0.4 rad, 
the two fundamental strings are splitted. In the future, it 
will be interesting to complement the present study of the 
flux tubes, with the computation of the static potential 
for the U geometry. 

ACKNOWLEDGMENTS 

This work was financed by the FCT contracts 
POCI/FP/81933/2007 and CERN/FP/83582/2008. We 
thank Orlando Oliveira for useful discussions and for 
sharing gauge field configurations.. 



[1] R. W. Haymaker, V. Singh, Y.-C. Peng, and J. Wosiek, [3] H. D. Trottier, Nucl. Phys. Proc. Suppl. 47, 286 (1996), 

Phys. Rev. D53, 389 (1996), arXiv:hep-lat/9406021. arXiv:hep-lat/9511006. 
[2] T. Barczyk, Acta Phys. Polon. B26, 1347 (1995). 



7 





(a) d = and 1 = 8 (b) d = 2 and Z = 8 (c) d = 4 and ■ 



(d) d = 6 and Z = 8 





-4 -2 2 4 



(e) n = 2 and r2 = 8 (f) n = 4 and r2 



(g) n. = 6 and r2 = 8 



(h) ri = 8 and r2 = 8 



Figure 9: Results for the lagrangian density. The Fig. (a)-(d) are for the U shape geometry and the Fig. (e)-(h) are 
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